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Abstract 

We evaluate the Casimir energy and force for a massive scalar field with general curvature 
coupling parameter, subject to Robin boundary conditions on two codimension-one parallel 
plates, located on a (D + l)-dimcnsional background spacetime with an arbitrary internal 
space. The most general case of different Robin coefficients on the two separate plates is 
considered. With independence of the geometry of the internal space, the Casimir forces are 
seen to be attractive for special cases of Dirichlet or Neumann boundary conditions on both 
plates and repulsive for Dirichlet boundary conditions on one plate and Neumann boundary 
conditions on the other. For Robin boundary conditions, the Casimir forces can be either 
attractive or repulsive, depending on the Robin coefficients and the separation between the 
plates, what is actually remarkable and useful. Indeed, we demonstrate the existence of an 
equilibrium point for the interplate distance, which is stabilized due to the Casimir force, 
and show that stability is enhanced by the presence of the extra dimensions. Applications of 
these properties in braneworld models are discussed. Finally, the corresponding results are 
generalized to the geometry of a piston of arbitrary cross section. 



PACS numbers: 03.70. +k, ll.10.Kk 



1 Introduction 

Many of the high-energy theories of fundamental physics are formulated in higher-dimensional 
spacetimes. In particular, the idea of extra dimensions has been extensively used in supergravity 

*Also at Center of Theor. Phys., TSPU, Tomsk 
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and superstring theories. It is commonly assumed that the extra dimensions are compactified. 
From the inflationary point of view, universes with compact spatial dimensions, under certain 
conditions, should be considered a rule rather than an exception [1]. Models involving a com- 
pact universe with non-trivial topology play a very important role by providing proper initial 
conditions for inflation. And compactification of spatial dimensions leads to a number of inter- 
esting quantum field theoretical effects, which include instabilities in interacting field theories, 
topological mass generation, and symmetry breaking. 

In the case of non-trivial topology, the boundary conditions imposed on fields give rise to a 
modification of the spectrum for vacuum fluctuations and, as a result, to Casimir-type contribu- 
tions in the vacuum expectation values of physical observables (for the topological Casimir effect 
and its role in cosmology see [2] and references therein). In models of the Kaluza-Klein type, 
the Casimir effect has been used as a stabilization mechanism for moduli fields and as a source 
for dynamical compactification of the extra dimensions, in particular, for quantum Kaluza-Klein 
gravity (see Ref. [3]). The Casimir energy can also serve as a model for dark energy needed 
for the explanation of the present accelerated expansion of the universe (see [3] and references 
therein). In addition, recent measurements of the Casimir forces between macroscopic bodies 
provide a sensitive test for constraining the parameters of long-range interactions, as predicted 
by modern unification theories of fundamental interactions [5] . The influence of extra compact- 
ified dimensions on the Casimir effect in the classical configuration of two parallel plates has 
been recently discussed in [6]- [9], for the case of a massless scalar field with Dirichlet bound- 
ary conditions, and in |10j-|13j. for the electromagnetic field for perfectly conducting boundary 
conditions. 

More recently, interest has concentrated on the topic of the Casimir effect in braneworld 
models with large extra dimensions. In this type of models (for a review see [13]) the concept 
of brane is used as a submanifold embedded in a higher dimensional spacetime, on which the 
standard-model particles are confined. Braneworlds naturally appear in the string/M theory con- 
text and provide a novel set up for discussing phenomenological and cosmological issues related 
with extra dimensions. In braneworld models the investigation of quantum effects is of consid- 
erable phenomenological interest, both in particle physics and in cosmology. The braneworld 
corresponds to a manifold with boundaries. All fields which propagate in the bulk will give 
Casimir-type contributions to the vacuum energy and, as a result, to the vacuum forces act- 
ing on the branes. Casimir forces provide a natural mechanism for stabilizing the radion field 
in the Randall- Sundrum model, as required for a complete solution of the hierarchy problem. 
In addition, the Casimir energy gives a contribution to both the brane and the bulk cosmo- 
logical constants. Hence, it has to be taken into account in any self-consistent formulation of 
the braneworld dynamics. The Casimir energy and corresponding Casimir forces within the 
framework of the Randall-Sundrum braneworld [15] have been evaluated in Refs. |16j-|22| by 
using both dimensional and zeta function regularization methods. Local Casimir densities were 
considered in Refs. \23\ 124] . The Casimir effect in higher dimensional generalizations of the 
Randall-Sundrum model with compact internal spaces has been investigated in [25j-[30j. 

The purpose of the present paper is to study the Casimir energy and force for a massive scalar 
field with an arbitrary curvature coupling parameter, obeying Robin boundary conditions on two 
codimension one parallel plates which are embedded in the background spacetime R( Dl ~ 1 > 1 > x S, 
being £ an arbitrary compact internal space. The most general case is considered, where the 
constants in the boundary conditions are different for the two separate plates. It will be shown 
that Robin boundary conditions with different coefficients are necessary to obtain repulsive 
Casimir forces. Robin type conditions are an extension of Dirichlet and Neumann boundary 
conditions and genuinely appear in a variety of situations, including vacuum effects for a confined 
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charged scalar field in external fields [31], spinor and gauge field theories, quantum gravity and 
supergravity [32j. Robin conditions can be made conformally invariant, while purely-Neumann 
conditions cannot. Therefore, Robin type conditions are needed when one deals with conformally 
invariant theories in the presence of boundaries and wishes to preserve this invariance. It is 
interesting to note that a quantum scalar field satisfying Robin conditions on the boundary 
of a cavity violates the Bekenstein's entropy-to-energy bound near certain points in the space 
of the parameter defining the boundary conditions 33J. Robin boundary conditions are an 
extension of those imposed on perfectly conducting boundaries and may, in some geometries, 
be useful for modelling the finite penetration of the field through the boundary, the skin-depth 
parameter being related to the Robin coefficient |341l35j . In other words, those are the boundary 
conditions which are more suitable to describe physically realistic situations. This type of 
boundary conditions naturally arise for scalar and fermion bulk fields in the Randall-Sundrum 
model [181 EH ES] and the corresponding Robin coefficients are related to the curvature scale and 
to the boundary mass terms of the field. Robin boundary conditions also appear in the study 
of Casimir forces between the boundary planes of films (for a recent discussion with references 
see, for instance, [37J). The Casimir effect in the geometry of two parallel plates with Robin 
boundary condition was investigated in Refs. [351 [381 EH HOI SI] • Note moreover that boundary 
problems with non-local boundary conditions can also be reduced to corresponding ones with 
Robin conditions, with the coefficients depending on the wave vector components along the 
plates [42]. 

The outline of the paper is as follows. In the next section we will consider the geometry of 
the problem and the corresponding eigenfunctions. The Casimir energy for two parallel plates in 
the general case for the internal subspace is evaluated in Sect. El The boundary-free and single 
plate parts will be extracted in a cutoff independent way. Applications to braneworlds are then 
discussed. In Sect.S]we consider the Casimir forces and show that depending on the coefficients 
in the boundary conditions these forces can be either attractive or repulsive. The asymptotic 
behavior of the forces, for small and large interplate distances, is given. As an application of the 
general results, in Sect. [5] a simple example with an internal space S 1 is discussed, for general 
periodicity condition along the compactified dimension. For this special example, we also present 
the boundary-free part and extract from the single plate parts the topological contributions. The 
corresponding generalizations for the internal spaces (S 1 ) 1 ^ and S are also given in detail. In 
Sect. El we extend the results for the Casimir energy and force to the case of the geometry of a 
piston with arbitrary cross-section. Section [7] contains a summary of the work. 

2 Geometry of the problem and eigenfunctions 

We consider a scalar field <p(x), with arbitrary curvature coupling parameter satisfying the 
equation of motion 

G? MJV V M Vjv + m 2 + (R) ip(x) = 0, (1) 

M, N = 0, 1, . . . , D, with R being the scalar curvature for a (D + l)-dimensional background 
spacetime (for the metric signature and the curvature tensor we adopt the conventions of Ref. 
|43] ) . For the special cases of minimally and of conformally coupled scalars one has, respectively, 
C = and £ = Cd = (D — 1)/4D. We will assume that the background spacetime has a topology 
^>(Di,i) x yj^ w j iere ft(Di,i) j g ^jj^ _|_ i )-<Timensional Minkowski spacetime and E a Z^-dimensional 
internal manifold, D = D\ + Di- The corresponding line element has the form 

ds 2 = g MN dx M dx N = 7]^dx^dx u - yud^dX 1 , (2) 



3 



with r)n V = diag(l, — 1, . . . , — 1) being the metric for the {D\ + l)-dimensional Minkowski space- 
time and the coordinates X" 1 cover the manifold S. Here and below fj,, v = 0,1,..., D\ and 
i, I = 1, . . . ,Z?2- For the scalar curvature of the metric tensor, from ^fy one has R = —R^, 
where Rr y ) is the scalar curvature for the metric tensor 7^. 

Our main interests in this paper will be to study the Casimir energy density and the mutual 
forces occurring for the geometry of two parallel infinite plates of codimension one, located at 
and x Dl = 02, ai < a2- As most general set up, we assume that on these boundaries 
the scalar field obeys Robin boundary conditions 

(1 + (3 jn M V M ) y(x) = [1 + fyi-iy^do.Mx) = 0, x Dl = a j} j = 1, 2, (3) 

with constant coefficients j3j. For (3j = these boundary conditions are reduced to Dirichlet one 
and for (3j = 00 to Neumann boundary conditions. The choice of different boundary conditions 
on the plates may correspond physically to use of different materials for plates. The imposition 
of boundary conditions on the quantum field changes the spectrum for the zero-point fluctua- 
tions and leads to the modification of the vacuum expectation values for physical quantities, as 
compared with the same situation without boundaries. 

In the region between the plates, a\ < x Dl < a>2, the corresponding eigenfunctions, satisfying 
the boundary condition on the plate at x Dl = aj, can be expressed in the decomposed form: 

I \ 

ip a (x M ) = C a exp -i V^x u cos [k Dl \x Dl - aj \ + a,] ipp{X), (4) 

\ M,^=0 / 

where a denotes a set of quantum numbers specifying the solution and 



k° = uj = ^k 2 + (k D i) 2 + m 2 p, mJ = Aj + m 2 , 

k = |k|, k= (fe 1 ,...,^ 1-1 ). (5) 



In Eq. ([3]), the ay, j = 1,2, are defined by the relations 

1 k^fa 



(6) 



The modes ipp(X) are the eigenfunctions of the operator A( 7 ) + C-R( 7 ) : 

[a (7) + CR(-y)] M x ) = ->$MX)> ( ? ) 

a, and fulfil 

r d D *X ^ p {X)iP},{X) = Spy. (8) 



with eigenvalues A|, and fulfill the normalization condition 



In Eq. d?]), A( 7 ) is the Laplace-Beltrami operator for the metric ju. In the consideration below 
we will assume that \p ^ 0. 

From the boundary condition on the second plate one obtains that the eigenvalues for k Dl 
are solutions of the equation 

F{z) = (1 — bib2Z 2 ) sin z — (61 + 62)2 cos z = 0, 

z = ak Dl , a = a,2 — 01, bj = Pj/a. (9) 
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We denote by z = z n , n = 1,2,..., the zeros of the function F{z) in the right half-plane of the 
complex variable z, arranged in ascending order, z n < z n +\. In the discussion below we will 
assume that all these zeros are real. This is the case for the conditions (see [3H]) {b\ + 62 ^ 
1,6162 ^ 0} U {6i 5 2 S; 0}. The coefficient C a in ([!]) is determined from the orthonormality 
condition for the eigenfunctions, and is equal to 

(10) 



being u(z n ) = Jk 2 + z 2 /a 2 + m 2 g the eigenfrequencies. 



(2VT) 1 

L0(z ri 



1 H sm(z n ) cos(z n + 2aj 

Zn. 



3 The Casimir energy 

The vacuum energy in the region between the plates (per unit volume along the directions 
x 1 , . . . , is given by the formal expression 

E M = \ J E E ^ 2 + z 2 ,a 2 + m 2 . (11) 

In the discussion below we will assume that some cutoff function is present, without writing it 
explicitly. Alternatively, one can use zeta function regularization, that yields the same result. 
For the sum over n we use the summation formula [381 IHj 

y ^f(Zn) = 7T /(0) | f°° 

1 + s'm(z n ) cos(z n + 2atj)/z n 2 1 - 6 2 - h J Q 



, f(iz) — f(—iz) , , 

+* / dz ( h J, u — ■ 12 

/ (biz-l)(b 2 z-l) 2z _ 1 v ' 

(6i^+l)(6 2 z+l) 



By taking into account the relation 



1 + ^ cos(z„ + 2a 3 ) = 1 - E -rA-i, (13) 

we see that the sum on the left-hand side of (|12|) coincides with the corresponding sum in the 
vacuum energy, if we take 

6,- 



/(*) = v fc2 + z2 ' a2 + m ly - E r^p] • (14) 

The application of the summation formula (|12p with (|14p allows us to write the vacuum 
energy from (jlip in the decomposed form 

^[01,02] = ^RPi-iJxE + E ^ + ^[01,03]' (I 5 ) 

3=1,2 

where we have introduced the notations 



^ = -57 p^?V fc2 + m ?-^/(Sp^i:y o d * r+ag ■ (16) 



/3 " " 7 
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and 



ai,a 2 \ 



AEi 



In Eq. flS 



7T 



dk 



(27T) 



/3 



/c 2 



(ftz+lXftz+l)' 



,2clz 



+E 



ft 



D 



(2») 



:E 

/3 



m. 



(17) 



(18) 



is the vacuum energy (per unit volume along the directions x 1 , . . . ,x Dl ) in the spacetime of 
topology i?^ 1 ' 1 ) x S for the case when the plates are absent. In the limit a — > oo the term 
A£ , [ 0l)Cl2 ] vanishes and the contribution Ej can be interpreted as the vacuum energy (per unit 
volume along the directions x 



induced by the presence of the plate located at 
dj in the half-space x Dl ^ aj. These single plate components do not depend on the 
location of the plate and do not contribute to the vacuum force acting on the plates. As it will 
be shown below, the latter is determined by the term Ai£r oa i. Note that this contribution is 
finite and that the cutoff function is strictly necessary for the terms E R (d 1 ,i) x ^ and Ej only. In 
the discussion below we will refer to A-Er ai)02 i as the interaction term. 

For further simplification of the corresponding expression, we use the relation 

dk r°° , r~, ~ ~ , , (4vr) 1 - Dl / 2 



(27T) 



dz 



k 2 



m 



2D 1 T(D 1 /2) 



dx(x 2 -m|) Dl / 2 5 (x). (19) 



In order to derive this formula, we must first integrate the left-hand side over the angular part 

I. After 



k 2 



of the vector k and then change to a new integration variable, y = w z" — — m^. 

introducing polar coordinates in the (k, y)-plane and integrating over the polar angle, we get 
Eq. (|19p . By using this relation, for the interaction part of the vacuum energy we find 



AEi 



(4tt) 



-£>i/2 



01,02 



r(Di/2 + i) 



E 



dx 



m 



!) Di/2 



(J3ix+l)(fax+l) 



!) e 2ax 



Pi 



(3 x 2 



(20) 



Using 



i 



(f3ix-l)(/3 2 x-l) Pax 
(hx+l^fox+l^ 



1 



1* 

dx 



1 



(fax - l)((3 2 x - 1) c _ 2ax 
(fax + l)((3 2 x + l) 



(21) 



and integrating by parts, the interaction term in the vacuum energy can be written as 

(4vr)- D i/ 2 



AE,, 



V{D x /2) 



/ dxx(x 2 



m 



2\ Dj/2-1 
/3J 



x In 



1 



(/?is + l)(fos + l) 2aa 
(/3ix-l)(/3 2 x-l) 



(22) 



In the case when the internal space is absent and for a massless scalar field this result reduces to 
the one derived in [35] . Note that the bulk divergences in the vacuum energy between the plates 
are contained in the first term on the right-hand side of (|15p and the boundary divergences are 
contained in the single plate contributions Ej. The interaction part is unambiguously defined. 
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In particular, it does not depend on the regularization scheme used (see, for example, Ref. 
for the case without the internal space, where exactly the same result is obtained with zeta 
function techniques) . 

For the special cases of Dirichlet and Neumann boundary conditions on both plates, from 
(p2l) one finds 

AE <"» , = _ <£ )- Di/2 ° y r d j* 2 -^ D "\ (23) 



with J = D and J = N for Dirichlet and Neumann boundary conditions, respectively. By ex- 
panding the factor l/(e 2ax — 1) in the integrand one gets 

,00 ( x 2_ m 2)D 1 /2 r(D/ 2 + l) " /flm ^(D 1+ l)/2 

e^-l = ^ g^) K (DM)/2 (2nam ), (24) 

being K u (z) the Mac-Donald (or modified Bessel) function. This allows us to write the corre- 
sponding vacuum energy for Dirichlet and Neumann scalars as 

Ap (J,J) 2a ~ Dl f(D 1+ i)/2^nam p ) 

^[ai.Oa] (Stt) (Di +l)/2 Z^Z^ n Di+l ' [ 0) 

y ' j3 n=l 

with the notation 

f v {z) = z v K v {z). (26) 

The energy given by Eq. (|25p is always negative and the corresponding Casimir forces are 
attractive for all interplate distances, as will be shown below. In the case D\ = 3 and for a 
massless scalar field, Eq. (|25|) reduces to the expression given in Ref. [SJ, where the zeta function 
method was used. 

For Dirichlet boundary conditions on one plate and Neumann boundary conditions on the 
other, similarly to (|20p we get 

AE (D,N) = W-*'*a V r (x*-m*)W 
rpi/2 + 1) yj mp e^ + l 

_ 2q- gl ^ ^ f{D 1+ \)/2^namp) 

( 8vr ) (Di +l)/2 Z^Z> (_l)n n £»i+l ' ^ ^ 

v ' /3 n=l v ' 

In this case the energy A£7r ai)02 i is always positive and the corresponding vacuum forces are 
repulsive for all distances between the plates. 

By using the result (|22p . in a similar way as in [21], we obtain the corresponding Casimir 
energy for a conformally coupled massless scalar field Tp(x) on the background of a spacetime 
with metric tensor g^fpf = Q 2 (x Dl )gMN, where the metric guN is defined by the line element 
d2|). We assume that the field obeys the boundary conditions: 

(I+^./Vm)*) = [1 + {-l)i- x toJ%d Dl ]Tj>{x) = 0, (lj = n(xf>), (28) 

on two codimension-one branes with coordinates x Dl = a,-, j = 1,2. The corresponding results 
for the interaction part of the Casimir energy can be derived from those obtained before simply 
by using the conformal relation that relates the two problems. The fields are connected by the 
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formula <p(x) = f^ 1 D ^ 2 (p(x). Making use of this relation, from Eqs. @ and (|28p we obtain the 
following relations between the Robin coefficients: 

-l 

Pj, (29) 

where Q'j = Vt'^x 1 ^ 1 ). We conclude that for a conformally coupled massless scalar field with 
boundary conditions (|28[) . the interaction part of the vacuum energy in the region between the 
branes is given by (|22j) . where the coefficients (5j are defined by the relations (|29p . In particular, 
for the case of Neumann boundary conditions = 0), one has (3j = 2Q,j{—iy/[{D — \)Q,'-]. In 

the special case of the AdS bulk used in the Randall-Sundrum braneworld model [15] (note that 
in this model only the inside region between the branes is considered) we have £l(x Dl ) = rr>/x Dl , 
being ro the AdS radius. The corresponding Robin coefficients for an untwisted scalar are given 
by the relations [H EH ES [36] 

pJ 1 = (-l) j c j /2-2D(/r D , (30) 

where c\ and C2 are the mass parameters in the surface action of the scalar field for the left and 
right branes, respectively. For a twisted scalar field, Dirichlet boundary conditions are obtained 
on both branes. 

To summarize, as we see, in the case of the warped geometry the corresponding vacuum 
energy is not, in general, a monotonic function of the inter-brane distance and can display a min- 
imum, corresponding to the stable equilibrium point. This property can be used in braneworld 
models for the stabilization of the radion field. An important difference between the warped 
geometry and the one discussed before is that now the single brane contributions to the vacuum 
energy depend on the location of the brane and, hence, give additional contributions to the force 
acting on the brane. The divergences in the single brane components are absorbed by adding 
the respective counterterms to the brane action. The coefficients of these counterterms are not 
computable within the framework of the low-energy effective theory and should be considered 
as parameters which are fixed by imposing renormalization conditions on the corresponding 
effective potential (see also the discussions in Refs. [TBI [IK] 1151 [251 I26] 1251 145]). 



Pi 



n 



4 The Casimir force 

The vacuum energy corresponding to the region ^ x l ^ q, I = 1, . . . , D\ — 1, a\ ^ x Dl a 2 
will be denoted £'r ai)a2 ici • • • cd^-i- The volume of this region is V = Vy,c\ ■ ■ ■ cd 1 -icl, being Vs 
the volume of the internal space. The corresponding vacuum stress at x Dl = a\+ is given by 

P(ai+) = —^E [aiMCl ■ ■ ■ c Dl _i = P + AP(ai+), (31) 

where 



Po = ~ , AP(a 1+ ) = - — ^E [aiM . (32) 



Using Eq. ([22]), we find 



AP( ai +) = -=^l -V dxx 2 -r-— jjr-JPL . (33) 

p 11 (p\x J rl)(fiix- s r\) 
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As can be easily seen, the vacuum stress at x Dl = a>2— is given by the same expression: 
AP(oa-) = AP(oi+). 

For the geometry of two parallel plates, the total vacuum energy is the sum of the contribu- 
tions from the regions x Dl ^ a%, a\ x Dl ^ 02 and ci2 ^ x Dl . When investigating the resulting 
force on the plate at x Dl = ax, in order to deal with finite regions from both sides, we will 
consider a piston like geometry (with large transverse dimensions, for a piston with finite cross 
section see below) , assuming the presence of an additional plate located at x Dl = ao < ax • For 
the corresponding vacuum stress at x Dl = ax — , one has 

P(ax-) = P + AP( ai -), AP( ai -) = -^-^AE [Kai] , (34) 

with b = ax — ao. The resulting pressure on the plate at x Dl = ax is given by the difference 

P( ai ) = AP{ax+)- AP(ax-). (35) 

As we see, the contributions to the vacuum force coming from the term P$ are the same from the 
left and from the right sides of the plate, so that there is no netto contribution to the effective 
force. In the limit ao — * —00 one has AP(ax—) — * and the Casimir force acting on the plate 
at x Dl = aj, j = 1, 2, in the original two-plate geometry is given by the expression 

p _ 2(4vr)-^/ 2 v [°° * 2 (* 2 ~ a 2 m}) D ^ 

V^(Dx/2)aM J (M-iK^-i) ^ _ 1 ■ W 

p p (b±x+X)(b2X+X) 

This force is attractive when P < and repulsive when P > 0. If one does not take into account 
the contributions from the exterior regions x Dl ^ ax and x Dl ^ 02 , the effective pressure is given 
by (|3ip where the renormalized value for Pq does not depend on the separation of the plates. 
In the case £ = S 1 and for periodic boundary conditions along the compactified dimension, the 
renormalized value Pq is positive (see next section, Eq. (I47p ) which would correspond to the 
repulsive force between the plates observed in the first paper of j6] (see also the discussion in 
Ref. 0). 

As is clearly seen from Eq. (j36f) . the sign of the vacuum stress on the plate is determined 
by the sign of the integral in this formula. In Fig. Q] we have plotted the location of the zeros 
for this integral on the (61, f)2)-plane in the case Dx = 3 and for different values of amp (figures 
on the curves). For a given amp, these zeros are located on two curves symmetric with respect 
to the line bx =62- The integral is positive in the region containing this line and it is negative 
outside. In particular, the Casimir force between the plates is always attractive for symmetric 
boundary conditions with fix = P2 < 0. This result is a special case of the general theorem [46], 
which dictates an attraction between bodies with the same properties. Note that the curves in 
Fig. [T] display the locations of the zeros for the Casimir force in the geometry of two parallel 
plates on the background of a 4-dimensional Minkowski spacetime. 

For the special cases of Dirichlet and Neumann boundary conditions, making use of the 
recurrence relations for the function K u (z), the Casimir forces can be written as 

= 2(4,)-^ v r d xx ^ x2 - m ^ )Dl12 " 

V s T(Dx/2) ^ J m0 e*o*-l 

EE rol+l [f{D 1 +x)/2(2namp) f(D 1 +3)/2(2namp)] , (37) 
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Figure 1: The location of the zeros for the integral in 
values of the parameter amp (numbers near the curves). 



in the case D\ = 3 and for different 



with J = D, N. Again, in the case D\ = 3 and for a massless scalar field this result reduces to the 
one obtained in [8]. The forces described by Eq. (|37p are attractive for all distances between the 
plates, irrespective of the geometry of the internal subspace. For Dirichlet boundary conditions 
on one plate and Neumann boundary conditions on the other, the expression for the Casimir force 
is obtained from Eq. (|37p after introducing an additional factor (— l) n+1 in the summation over 
n. In this case, Casimir forces are repulsive. In the general case of Robin boundary conditions 
the Casimir forces can be either attractive or repulsive, depending on the coefficients present 
in the definition of the boundary conditions, and on the distance between the plates. For the 
special case of the topology R^ D ~ 1 ^ x S 1 this issue will be illustrated in the next section. 

Let us now consider the asymptotic behavior of the Casimir force as a function of the size of 
the internal space. Note that if the size of the internal space is of the order L, then for nonzero 
modes one has \p ~ 1/L. For small values of L and for the nonzero modes, A^ is large. The 
contribution of these modes is exponentially suppressed and the main contribution comes form 
the zero mode. In this case, from (|33|) we recover the Casimir force for two parallel plates in 
(D\ + l)-dimensional Minkowskian spacetime, namely 



2(4vr) 



-Dl/2 



rpi/2) 



dx 



x 2 (x 2 



m 



!\ d/2-l 



(.Pix-l)(lhx-l) 2ax 



1 



(38) 



For the case of a degenerated zero eigenstate the corresponding degeneracy factor must be 
included on the right-hand side. In some models of compactification the zero mode is absent 
(for example, in models with twisted boundary conditions along the compactified dimensions, 
see below). In such cases, for small values of L the main contribution to the Casimir force comes 
from the lowest mode \p = Aq and, to leading order, one gets 



in 



£>i+l 




(Pim x + l)(/?2"i x + 1) e 



-2amo 



(4tt) d i/ 2 (Pimox - l)(p 2 m x - 1) (am ) D ^ 2 ' 



(39) 



where mo = \/Ag + ni 2 . Hence, here the Casimir forces are exponentially suppressed for small 
sizes of the internal space. 
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For small values of the inter-plate distance, a/\/3j\ <C 1, the main contribution into the 
integral in Eq. f)33|) comes from larger values of x and, to leading order, one has 

- T ; =f= — > / dxx 2 - , 40 

F E r(L» 1 /2)^J m ^ e 2 --l V ; 

except for the case of Dirichlet boundary conditions on one plate and non-Dirichlet boundary 
conditions on the other. We see that in this limit (|40p the Casimir force is attractive. However, 
in the case of Dirichlet boundary condition on one plate and non-Dirichlet boundary condition 
on the other the Casimir force is repulsive at small distances, what is indeed a remarkable result. 

5 Particular cases 

As a simple example of a particular application of the general results obtained above, we will 
first consider the special case where £ = S , with the size of the internal space being 2-kL. For 
the compact dimension we assume a general periodicity condition of the form 

MX + 2ttL) = e 2 ™V /3 (X), (41) 

with constant a, ^ a ^ 1. The specific cases a = and a = 1/2 correspond to untwisted and 
to twisted fields, respectively. The corresponding part of the scalar eigenfunctions is 

iKX 

MX) = ^==, K = (P + a)/L, /3 = 0, ±1, ±2, . . . . (42) 

The formulas for the topology i^- 0-1 ' 1 ) x S 1 are obtained from the results given in the previous 
sections, by taking 

+ OO 1/3 -L. I 

E = E > A /3 = —j^-, mp = y/^ + a) 2 /L^ + m 2 , D x =D-l. (43) 

P P=-oo 

In particular, for the Casimir force one has 

(Air)-^- 1 ^ 2 +^ /•«> x 2 (x 2 -m2)(D-3)/2 
P = ~TrT((D-l)/2)L E J m dx W^)W 2 x-i) 2ax 7/ ^ 

For large values of L the main contribution to the series in (|44p comes from large values of [3 
and one can replace the summation over (3 by an integration. After some transformations, to 
leading order we find the result 

^~ rp/2) i ^-a^i e^-l' 

(/3ix+l)(/32a;+l) 

which is, in fact, the Casimir force for two parallel plates in (D + l)-dimensional Minkowskian 
spacetime. 

As we already noted before, depending on the values of the coefficients (3j and of the distance 
between the plates, the Casimir force (|44|) can be either attractive or repulsive. In Fig. [51 and 
corresponding to the model with D = 4 and for a massless scalar field with Dirichlet boundary 
conditions, we have plotted the ratio 2-kLP/Pc as a function of a/L, where Pq = — vr 2 /(480a 4 ) 
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1/2 



0.0 



a/L 



Figure 2: Ratio of the Casimir force for two parallel plates in the spacetime with topology 
Ri 3 * 1 ) x S 1 to the standard Casimir force in i^ 3 ' 1 ), for a massless Dirichlet scalar, as a function 
of a/L. The values on each of the curves correspond to those of the parameter a. 

is the standard Casimir force. The values on each of the curves correspond to those of the 
parameter a. As we have explained before, for a / the zero mode is absent and for large 
values of a/L the Casimir force is exponentially suppressed. 

In Fig. [3] the Casimir force is plotted for the topology i?( 3,1 ) x S 1 and an untwisted (a = 0) 
massless scalar field with Robin coefficients (3\/ao = —0.1, fo/ao = —0.5 (a a fixed length scale) 
as a function of a/a$ for L/a$ = 1 (full curve) and L/ao = 0.5 (dashed curve). The thick curve 
corresponds to the Casimir force for two parallel plates in Minkowski spacetime it^ 3,1 ) with the 
same Robin coefficients. As is seen, the corresponding Casimir forces are attractive for small 
and large distances between the plates while they are repulsive for intermediate distances. There 
are two equilibrium points corresponding to the zeros of the function P. The leftmost point is 
unstable whereas the rightmost one is stable. Hence, in this case the Casimir force stabilizes the 
distance between the plates. This feature can be used in braneworld models for the stabilization 
of the radion field. We see from Fig. [3] that the height of the barrier between the stable and the 
unstable equilibrium points is increased by the presence of the internal space. As a consequence, 
an enhancement of the repulsive Casimir effect, coming from the extra dimension, occurs. 

As already explained before, in the case of Dirichlet boundary conditions on one plate and 
non-Dirichlet boundary conditions on the other, the Casimir force is repulsive at small separa- 
tions. In Fig. H] we illustrate this feature for the topology I?' 3 ' 1 ) x S 1 in the case of an untwisted 
massless scalar field with /3i = 0, fo/a-o = —0.5. As in Fig. O the full (dashed) curve stands 
for L/ao = 1 (L/ao = 0.5) and the thick curve corresponds to the Casimir force for two parallel 
plates in Minkowski spacetime i?( 3,1 ) with the same Robin coefficients. 

For the topology under consideration, the Casimir energy for the bulk without boundaries, 

E R (d-i,i )xS i = ^J ^D-i + (/? + «) 2 /£ 2 + m2 
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Figure 3: Casimir force for two parallel plates in the spacetime with topology R^' 1 ^ x S 1 , for 
an untwisted massless scalar field with Robin coefficients fii/ao = —0.1, ^/oq = —0.5, as a 
function of the distance between the plates. The full (dashed) curve corresponds to a size of the 
internal space with L/ciq = 1 (L/oq = 0.5). The thick curve corresponds to the Casimir force 
for two parallel plates in Minkowski spacetime i?( 3,1 ) with the same Robin coefficients. 
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can be further simplified through the Abel-Plana summation formula, in the form [444 |4"7] 



E f(\P + a\)=2r dxf ( X ) +i r dx £ ^4^ , (46) 



/3=-oo ,/u - 70 A=±l 



what leads to the result 

2vrL /■ d D k 



E R( o-i,i )xS , = — 7^D^l+m 2 



2(Lm)( 0+1 V 2 ~ cos(27rna) T ^ tn r , 

' (27r d l d E n (D + i)/2 ^ (D + i)/ 2 (2vrnLm). (47) 

n=l 



The second term on the right-hand side of this expression is finite and introduction of a cutoff 
function is necessary for the first term only. Note that the latter is the vacuum energy density 
for the spatial topology R D and, hence, the second term on the right-hand side of Eq. (|4Tj) is 
the contribution to the vacuum energy induced by the compactness of the x D dimension. In 
particular, the topological part of the vacuum energy is always negative (positive) for untwisted 
(twisted) scalars. 

In a similar way, we can also extract the topological contributions in the single plate terms 
of the vacuum energy. After applying the summation formula (|46p and after integration, these 
contributions yield 

n=l 

(4ir)-(*-i>/2Lft v r» (y 2 -mT /2 f 1 , (1 - s 2 )^/ 2 ^ 

r((D + l)/2) Z^J m aV ^L V+ 2niXa_ 1 J 1 + ^2 (y2 _ m 2 )x 2 > ^ 

where 



(27r)0 -iV^-l^"' 27r7 (2vr)«-iy "* 1 + /3|x 2 

is the vacuum energy (per unit volume along the coordinates x 1 , . . . , re 0-1 ) for a single plate in 
Minkowski spacetime with trivial topology R^ D,1 \ Hence, the last two terms on the right-hand 
side of (|48p are the terms in the vacuum energy corresponding to a single plate and due to the 
compactness of the dimension x . These terms are finite and renormalization is needed for the 
Minowskian part Ej^ only. Note that for Dirichlet and Neumann boundary conditions the last 
term on the right of Eq. (|48p vanishes. 

The case of a /^-dimensional torus as internal space, E = (S l ) D ' 2 , can be considered in a 
similar way. For a scalar field with the periodicity condition ipp(X l + 2ttLi) = e 2,Kioil ipg(X ) 
along the coordinate X 1 , ^ X 1 ^ 2ttLi, the formulas for the Casimir energy and force are 
obtained from the general expressions in Sects. [3] and H] with the substitutions 

+oo +oo £>2 

E= E ■- E > ™} = EUi + <*?/lZ + m*. (50) 

fi ji=-oo jD 2 =-° i=l 

Concerning the issues of embedding the model in string theory and of the discussions of the 
holographic principle there, the case of the internal space S = S Di is of very special interest. 
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The corresponding eigenfunctions ijjg(X) are expressed in terms of spherical harmonics of degree 
I, I = 0, 1, 2, . . .. For the internal space with radius L the expressions for the Casimir energy and 
Casimir force are quite easily obtained from the general formulas given in the above sections, 
just by replacing 

T{1 + D 2 -l) 



£ - 5> + i*-i)- 

p i=o v l ' 

1 



A/? -> -y/l(l + D 2 -l)+CD 2 (D 2 -l). (51) 
Here the factor under the summation sign is the degeneracy of the angular mode with a given I. 



6 Generalized piston geometry 

In a way very much similar to the procedure described in the preceding sections, we are able 
to treat the more general case when a part of the dimensions x 1 , . . . , are still constrained 

by boundary conditions. This corresponds to considering a generalized piston geometry, a quite 
fashionable situation nowadays, in particular, in the quest for negative Casimir forces (for the 
investigation of the Casimir effect in a piston geometry see |48j and references therein). Here 
we have obtained those in the configurations above, but would like to see now the differences 
introduced in our results by the consideration of piston geometries. 

We will denote by d\ the number of unconstrained dimensions (coordinates x 1 , . . . , x dl ) and 
by 7?, with a collective index i, the eigenvalues of the Laplacian along the constrained directions: 

& Dl -i-d^ a {x M ) = -jf^ a (x M ). (52) 
The eigenfrequencies in the region between the plates are here given by 

w(*0 = \/k 2 dl +zl/a 2 + jf + m 2 p, (53) 

with the vacuum energy being 

*W 2 ] = \f TotEE y/^+^ 2 + lf + m% (54) 
J y ' i,/3 n=l 

After applying the summation formula (|12p to the sum over n, we can write the energy in the 
decomposed form 

J K ' i,/3 j=l,2 

Here the expressions for the terms Ej and A£ , [ aij(l2 ] are obtained from Eqs. ()16f) and (11T[) by the 
replacements 

/ (2^i-i ~* / ~" (56) 

^ -»■ D x d x + 1, mj -» % 2 + mjj. (57) 

/3 i,/3 
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These formulas are further simplified after integrating over the angular part of the vector k^. 
The corresponding expressions are obtained from the results of Sects. [3]and[U with the replace- 
ments (|57p . In particular, for the interaction part of the Casimir energy (per unit volume along 
the direction x 1 , . . . , x dl ) one has 

AE, i = 1 ] — V / dxxix 2 -7 2 -mD^-W 



x In 



_ (f3 lX + l)(f3 2 x + l) 2ax 
(P lX - l)((3 2 x - 1) 



(58) 



The expression for the Casimir pressure takes the form 



2 U 7T X-(d 1 +l)/2 roc x 2^2_ 7 2_ m 2)(<i 1 -l)/2 

P(a ' ^ ^ = " V^Vmck + l)/2) 2l\ n^ dx (ft«-i)(ft*-i) c ^ — ' (59) 

where V^ s is the volume of the piston cross section along the coordinates x dl+l , . . . ,x Dl ~ l . In 
particular, for Dirichlet and Neumann boundary conditions on the plates, we find 

2a -di-2 oo 

P(J ' J) (°) = (8T)jl /2 + l Us EE^[^/w( Z ) - ^/2+2 W] z=2na y^|> (60) 

J = D,N, the function /^(.z) being defined in Eq. (|26p . The corresponding force remains attrac- 
tive independently of the form of the cross section. In the case of Dirichlet boundary condition 
on one plate and Neumann boundary condition on the other the expression for the Casimir force 
is obtained from (|60f) after introducing an additional factor (— l) n , and the resulting force is 
always repulsive. In the special case d% = and for a massless scalar field, Eq. (|60p reduces to 
the formula for the Casimir force in Ref. [;8J, as it should. 

On the base of Eqs. ()58[) and (|59p we can analyze the geometry of a generalized piston with 
two chambers, assuming that the plates are located at x Dl = ao, ai, a 2 , with the Robin coefficient 
(3q for the left plate (for the Casimir effect in a piston geometry see, for example, Refs. |48|). 
For an arbitrary cross section, the effective pressure on the plate at x Dl = a\ is given by 

P ai =P(a,p 1 ,fo)-P(b,p ,p 1 ), (61) 

with b = a\ — ao. For P ai < (P ai > 0) the resulting force on the plate is directed towards the 
right (left) plate. For Dirichlet and Neumann boundary conditions the Casimir stress given by 
Eq. (|60p is a monotonic function of the plate separation and, hence, in the piston geometry with 
two chambers the resulting force (|6ip is directed toward the closer plate. 

In the special case of the geometry of a piston with circular cross section in the plane 
(x Dl ~ 2 , x^* 1-1 ), the corresponding part in the eigenfunctions has the form J\ q \{j]r)e iq ^ , q = 
0, ±1, ±2, where r and are polar coordinates on this plane and J q (z) is the Bessel function. 
The eigenvalues for the quantum number rj are quantized by the boundary conditions on the 
cylindrical surface r = ro, with vq being the piston radius. For example, in the case of a 
Dirichlet boundary condition one has rj = j\ q \ :P /ro, where j UtP is the p-th positive zero of the 
function J v {z). The corresponding formulas for the Casimir energy and forces are obtained from 
the general results (f58f) and (|36|) with the substitutions 

+oo oo 

di = Di - 3, = IE' 7* = j\q\, P / r o- ( 62 ) 

i q=—oop=l 
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Other types of boundary conditions on the cylindrical boundary can be considered in a similar 
way. Moreover, the generalization of our procedure to more than three plates is also easy to 
carry out. 

7 Conclusion 

We have investigated in this paper the influence of extra dimensions on the Casimir energy and 
on the Casimir force for a massive scalar field with an arbitrary curvature coupling parameter, 
in the usual geometry of two parallel plates. We have assumed that on the plates the field obeys 
Robin boundary conditions with, in general, different coefficients for the two different plates. The 
corresponding eigenfrequencies are expressed in terms of solutions of a transcendental equation 
0, thus they are known implicitly only. By applying the summation formula (|12p to the 
corresponding series in the mode-sum for the vacuum energy in the region between the plates, we 
have explicitly extracted, in a cut-off independent way, the boundary-free (topological) part and 
the contributions induced by the single plates (when the other plate is absent). The remaining 
interaction part is finite for all nonzero inter-plate distances and is cut-off independent. The 
surface divergences in the Casimir energy are contained in the single plate components only. 
But the latter do not depend on the location of the plate and do not contribute to the Casimir 
force. For an arbitrary internal space, the interaction part of the Casimir energy is given by 
Eq. (|22p . In the special cases of Dirichlet and Neumann boundary conditions on both plates 
this formula leads to the result f)25f) and the corresponding energy is always negative. For 
Dirichlet boundary conditions on one plate and Neumann boundary conditions on the other, 
the interaction component of the vacuum energy is given by Eq. (|27p . and it is positive for all 
values of the interplate distance. In the case of a conformally coupled massless field on the 
background of a spacetime conformally related to the one described by the line element ([2]) with 
the conformal factor Q 2 (x Dl ), the interaction part of the Casimir energy is given by Eq. (|22h . 
with the coefficients j3j being related to the specific coefficients of the Robin boundary conditions 
(I28j) and to the conformal factor by Eqs. (I29h . In the Randall- Sundrum two brane model with 
a compact internal space, the corresponding Robin coefficients are given by Eq. (|30p and the 
corresponding vacuum energy can have a minimum, corresponding to the stable equilibrium 
point. This feature is useful in braneworld models for the stabilization of the radion field. 

The interaction forces between the plates for the most general case of internal space have been 
considered in Sect. 01 In order to obtain the resulting force, the contributions from both sides 
of the plates must be taken into account. Then, the forces coming from the topological parts of 
the vacuum energy cancel out and only the interaction terms contribute to the Casimir force. 
In order to show this important fact explicitly, we have considered a piston-like geometry, by 
introducing a third plate. At the end of the calculation this plate is sent to infinity. The resulting 
Casimir force is given by Eq. (|36|) . With independence of the geometry of the internal space, 
the force is attractive for Dirichlet or Neumann boundary conditions on both plates (formula 
(|37p ) and it is repulsive for Dirichlet boundary conditions on one plate and Neumann boundary 
conditions on the other. In both cases the force is a monotonic function of the distance. For 
general Robin boundary conditions the Casimir force can be either attractive (corresponding 
to negative values of P) or repulsive (positive values of P), depending on the particular Robin 
coefficients and on the distance between the plates. For small values of the size of the internal 
space and in models where the zero modes along the internal space are present, the main 
contribution to the Casimir force comes from the zero modes and the contributions of the nonzero 
modes are exponentially suppressed. In this limit, to leading order we recover the standard result 
for the Casimir force between two plates in (D\ + l)-dimensional Minkowski spacetime. When 
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the zero mode is absent (for example, in the case of twisted boundary conditions along the 
compactified dimensions), the Casimir forces are exponentially suppressed in the limit of small 
size of the internal space. For small values of the inter-plate distance the Casimir forces are 
attractive, independently of the values of the Robin coefficients, except for the case of Dirichlet 
boundary conditions on one plate and non-Dirichlet boundary conditions on the other. In this 
latter case, the Casimir force is repulsive at small distances. It is interesting to remark that this 
property could be used in the proposal of a Casimir experiment with the purpose to carry out 
an explicit detailed observation of 'large' extra dimensions as allowed by some models of particle 
physics. 

As an illustration of the general results, in Sect. [5] we have considered a special model for the 
internal space £ = S , with the periodicity condition (|4ip along the compactified dimension. 
For the specific values a = 0,1/2 this condition corresponds to untwisted and twisted scalar 
fields, respectively. In Fig. [2 for Dirichlet boundary conditions we depicted the dependence of 
the Casimir force on the distance between the plates for different values of the parameter a. In 
Fig. [3] a plot of the Casimir force in the case of Robin boundary conditions on both plates as a 
function of the inter-plate distance has been provided. In the example considered, the Casimir 
force is attractive both for large and for small distances, while it is repulsive at intermediate 
distances. The Casimir force vanishes at two values of the inter-plate distance, which correspond 
to equilibrium points. The leftmost point is unstable and the rightmost one is locally stable. As 
shown in the plot, the stability of the rightmost equilibrium point is enhanced by the presence 
of the internal space. Formulas for the Casimir energy and force for the more general internal 
spaces (S 1 ) 02 (ZVtorus) and S D2 have been obtained from the general results of Sects. [3] and HI 

In the last section [H we have extended the results from the previous ones to the case of 
a piston geometry with finite cross section of arbitrary form along some subset of the dimen- 
sions. The corresponding expressions for the Casimir energy and Casimir force between the 
plates are given by Eqs. ([58]) and ([59]) . We have checked that the qualitative features described 
above remain basically unaltered. In particular, the possibility of a repulsive Casimir effect is 
again observed. In the special case of a piston geometry of circular cross section on the plane 
the corresponding formulas have been specified in (|62j) . 

The search for specific applications of this study to practical situations in braneworld models, 
nano-physics and particle physics will keep us busy for some time. 
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